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Abstract. In this work we consider local minimizers (in the topology of transport distances) of the 
interaction energy associated to a repulsive-attractive potential. We show how the dimensionality 
of the support of local minimizers is related to the repulsive strength of the potential at the origin. 

(N 

^ 1. Introduction 

1^ Given a Borel measurable function W : — t- (—00, +00] which is bounded from below, the 

interaction energy of the Borel probability measure fi is given by 

^ Ew[fA-=l [[ Wix-y)df,{x)di^{y). (1) 

• Our main goal will be to analyse the qualitative properties of local minimizers of the energy Ey\r 

in the set of Borel probability measures with the topology induced by transport distances. More 
^ specifically, we will show that the Hausdorff dimension of the support of local minimizers is directly 

G related to the behavior at the origin of AW. 

The interaction energy Ew arises in many contexts. In physical, biological, and material sciences 
^— I it is used to model particles or individuals effects on others via pairwise interactions. Given n 

^ particles located at Xi, . . . , X^ G , their discrete interaction energy is given by 



in 

EUXi.....Xn]:=—^Y.^^^^-^^)- (2) 

O 

^SJ Note that formally for a large number of particles, the discrete energy ([2j) is well approximated by 

the continuum energy ([T| where dn{x) is a general distribution of particles at location x G M^. In 
J> fact, the continuum energy ([T]) of the discrete distribution ^ Y17=i reduces to ([2]). 

IT^ In models arising in material sciences [151 [371 [311 [211 [35] , particles, nano-particles, or molecules 

^ self-assemble in a way to minimize energies similar to Analogously in applications to biological 

sciences \27 \ \26 \ [33l HI IT^ . individuals in a social aggregate (e.g., swarm, flock, school, or herd) 
self-organize in order to minimize similar type of energies. In these applications the potential W is 
typically repulsive in the short range so that particles/individuals do not collide, and attractive in 
the long range so that the particles/individuals gather to form a group or a structure. Therefore 
one is often led to consider radially symmetric interaction potentials of the form W{x) = w{\x\) 
where w : [0, +00) — >• (—00, +00] is decreasing on some interval [0, vq) and increasing on (ro, +00). 
The function w may or may not have a singularity at r = 0. We will refer to such potentials as 
being repulsive-attractive. Since w has a global minimum at tq, it is obvious that if we consider 
only two particles Xi and X2, in order to minimize £^^[Xi,X2], the two particles must be located 
at a distance tq from one another. Whereas the situation is simple with two particles, it becomes 
very complicated for large number of particles. Recent works [231 [33 El [351 [22l [32l [301 (13 (IBl 13] 
have shown that such repulsive-attractive potentials lead to the emergence of surprisingly rich 



geometric structures. The goal of the present paper is to understand how the dimensionahty of 
these structures depends on the singularity of AW at the origin. 

Let us describe the main results. Consider a repulsive-attractive potential W{x) = w{\x\). 
Typically the Laplacian of such potential will be negative in a neighborhood of the origin. We show 
that if 

AW(x) ^ asx^O (3) 

\x\P 

for some < (3 < N , then the support of local minimizers of Ey\r has Hausdorff dimension greater 
or equal to (3. The precise hypotheses needed on W for this result to be true, as well as the 
precise meaning of ([3]), can be found in the statement of Theorem [T] The exponent /? appearing in 
^ quantifies how repulsive the potential is at the origin. Therefore our result can be intuitively 
understood as follows: the more repulsive the potential is at the origin, the higher the dimension 
of local minimizers will be. 

Potentials satisfying Q have a singular Laplacian at and we refer to them as strongly repul- 
sive at the origin. The second main result is devoted to potentials which are mildly repulsive 
at the origin, that is potentials whose Laplacian does not blow up at the origin. To be more 
precise we show that if 

W{x) ~ —\x\" as X — )• for some a > 2 (4) 

then a local minimizer of the interaction energy cannot be concentrated on smooth manifolds of any 
dimension except 0-dimensional sets. The exact hypotheses on W, as well as the precise meaning 
of 0, can be found in Theorem [2| Note that this result suggests that local minimizers of the 
interaction energy of mildly repulsive potentials have zero Hausdorff dimension - however we are 
currently unable to prove this stronger result. 

Summarizing, in this paper we show that if the Laplacian of the potential behaves like — l/jxl'^ 
around the origin, with < (3 < N, then the dimension of minimizers is at least /5 and if the 
Laplacian does not blow up at the origin, then the dimension is zero, see the precise statement in 
Theorems [T] and [2j This is illustrated in the case of two dimensions ( = 2) in Table [l| where we 
show some local minimizers of E\y with interaction potentials of the form 

W{x) = -— + — a<7, (5) 
a 7 

so that W{x) ~ and AW{x) ~ -yL with /? = 2- aasx->0. 

a \x\' 

• Subfigure (a): a = 2.5 and 7 = 15. The support of the minimizer has zero Hausdorff 
dimension in agreement with Theorem [2} Actually, in this particular case it is supported 
just on 3 points. 

• Subfigure (b) and (c): we consider two examples where the potentials have the same be- 
havior at the origin, a = 1.5, but different attractive long range behavior (7 = 7 and 
2 respectively) . Theorem [l] shows that the Hausdorff dimension of the support must be 
greater or equal to /3 = 2 — a = 0.5. Indeed, the minimizer for the first example has a 
one-dimensional support on three curves whereas the minimizer for the second example has 
two-dimensional support. 

• Subfigure (d): q = 0.5 and 7 = 5. Theorem [T] proves that the Hausdorff dimension of the 
support must be greater than /3 = 2 — a = 1.5. The numerical simulation demonstrates 
that it has dimension two. 

In our extensive numerical experiments using gradient descent methods we never observed mini- 
mizers with a support that might be of non-integer Hausdorff dimension. 

In most of this paper, we will consider local minimizers for the topology induced by the transport 
distance doo (see section 2 for a definition of doo)- This topology is indeed the natural one to 
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Table 1. Local minimizers of the interaction energy E^^ for various potentials 
W{x). In these computations n = 10,000. When ISW does not blow-up at the 
origin (Case a) the Hausdorff dimension of the the support of minimizers is zero. 
When ~ — l/|a;|^ asx— )-0, 0</3<A'' (Cases b,c,d) the Hausdorff dimension 
of the the support of minimizers is greater or equal to (5. 



consider. In particular, gradient descent numerical methods based on particles typically lead to 
local minimizers for the doo— topology. Moreover the topology induced by doo is the finest topology 
among the ones induced by dp, 1 < p < oo (see section 2 for a definition of dp). As a consequence 
local minimizers in the dp-topology are automatically local minimizers in the doo-topology, and thus 
they are also covered by our study. In Section 5 we will discuss in more detail these questions. 

Let us finally mention that the gradient flow of the energy E\\r in the Wasserstein sense d2 
[la H [13] has been extensively studied in recent years plElElISIHlIIIlIIOlElIITlIISlEOlElE]. It 
leads to the nonlocal interaction equation 

^ + div(/iv) = , v = -VW * n (6) 

where n{t,x) = ^t{x) is the probability or mass density of particles at time t and at location 
X € M'^, and v{t, x) is the velocity of the particles. Stability properties of steady states for ([g]) with 
repulsive-attractive potentials have only been analyzed very recently. In pj we gave conditions 
for radial stability/instability of particular local minimizers. We should also mention that the one 
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dimensional case was analyzed in detail in [TTl [18]. Well-posedness theories for these repulsive- 
attractive potentials in various functional settings have been provided in [JU [H [HI [TTl [SI [3] • Stable 
steady states of ([6) under certain set of perturbations are expected to be local minimizers of the 
energy functional (!]) in a topology to be specified. Actually, this topology should determine the 
set of admissible perturbations. As already mentioned, the doo -stability is the one typically studied 
by performing equal mass particles simulations. 

Finally, we can now interpret our dimensionality result in terms of the nonlocal evolution equation 
Q. The heuristic idea behind the implication: Q with < P < N implies dimensionality larger 
than /3 of the support of local minimizers of Eyy', can be understood in terms of the divergence of 
the velocity field in ([6]). In fact, it is straightforward to check that the divergence of the velocity 
field generated by a uniform density localized over a smooth manifold of dimension k is +oo on 
the manifold if and only if k < (3 (this is equivalent to non-integrability of — AVF on manifolds of 
dimension k). Heuristically, if div v = —AW*fi associated to fi diverges on its support the density 
has a strong tendency to spread, the configuration is not stable and then is not a local minimizer. 
Therefore, we can reinterpret our result in Theorem 1 as follows: local minimizers of ([T]) have to 
be supported on manifolds where the divergence of their generated velocity field is not +oo. 

The plan of the paper is as follows. Section 2 will be devoted to the necessary background in 
optimal transport theory and notations. Strongly repulsive potentials are treated in Section 3 while 
mildly repulsive potentials are analyzed in Section 4. In Section 5, for the smaller subset of local 
minimizers in the (i2-topology, we show that we can use an Euler-Lagrange approach in the spirit 
of [4J to derive some properties of these minimizers. Extensive numerical tests as well as details of 
the algorithm used in order to minimize E^r are reported in Section 6. 



We denote by i3(M^) the family of Borel subsets of M^. Given a set A G i3(M^), its Lebesgue 
measure is denoted by |^|. We denote by M{R^) the set of (nonnegative) Borel measures on M 
and by 'P(M^) the set of Borel probability measures on M^. The support of ^ G A4{M.^), denoted 
by supp(/i), is the closed set defined by 



A measure p G A^(M^) is said to be a part of p if p{A) < fi{A) for all A G B{M.^) and it is not 
identically zero. This terminology is justified by the fact that if /o is a part of p, then p can be 
written p = p + v for some v G A^(M^) {v = p — p to he more precise). We will say that a 
probability measure p G V{MJ^) can be decomposed as a convex combination of po,pi G P(M^) if 
there exists < mo, mi < 1 with mo + mi = 1 such that p = m^pQ + mipi. 



Let us introduce some notation related to the interaction potential energy. We denote by B\y : 
V{R^) X V{R^) (0, +oo] the bilinear form defined by 



Obviously we have that = B\y[p,p]- Let us define the shortcut notation T]y[pi,p2] := 

-^'VkI/^i] — "^BiyIpi, P2] + E\y[p2] which will often occur in several computations. For notational 
simplicity, we will drop the subscript for E\y, Bw, and T]y in detailed proofs while kept in the 
main statements. 

Let us give a brief self-contained summary of the main concepts related to distances between 
measures in optimal transport theory, we refer to [20l [25] for further details. A probability 
measure vr on the product space R^ x is said to be a transference plan between p G V{R^) 



2. Preliminaries in Transport Distances 



supp(/x) := {x G : piB{x, e)) > for all e > 0} . 




(7) 



and 1/ G V{R^) if 



tt{AxR'^) = p{A) and 7r(M^ x A) 



u{A) 



(8) 
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for all A £ B{R^). If E T'(M^), then 

n(/i, u) := {it G V{R^ xR^):^ holds for all A G B{R^)} 

denotes the set of admissible transference plans between /i and i'. Informally, if vr € n(/u, z^) 
then dTT{x,y) measures the amount of mass transferred from location x to location y. With this 
interpretation in mind note that supj-^. y■)gs^pp(^•) |a; — y| represents the maximum distance that an 
infinitesimal element of mass from n is moved by the transference plan vr. We will work with the 
oo-Wasserstein distance doo between two probability measures //, u defined by 

doo{lJ',v)= inf sup \x-y\, (9) 

7ren(^t,i/) {a;,j;)gsupp(7r) 

which can take infinite values, but it is obviously finite for compactly supported measures. This 
distance induces a complete metric structure restricted to the set of probability measure with finite 
moments of all orders, Voo{^-^), as proven in |2Dj. 

We remind that for 1 < p < oo the distance dp between two measures // and v is defined by 



<(m,i^)= inf {jl |x - y|Pd7r(x,y) 

Note that dp{^, v) < oo for fijf € Vp{R^) the set of probability measures with finite moments of 
order p. Since dp{fi, v) is increasing as a function of 1 < p < oo, one can show that it converges 
to doo(/^, z^) as p — )• oo. Since the distances are ordered with respect to p, it is obvious that the 
topologies are also ordered. More precisely, open sets for dp are always open sets for doo, and thus, 
doo induces the finest topology among dp, 1 < p < oo. More properties of the distance doo can be 
seen in [25j . 

Given T : — > R^ measurable, we say that u is the push-forward of fi through T, u = T#fi, 
if v[A] := n[T^^{A)] for all measurable sets A C R^ , equivalently 



/ 



ip{x)dv{x)= / ip(T{x)) d^{x) 



for all if G Cb(M^). In case there is a map T : R^ — > R^ transporting fj, onto u, i.e. T#/i = ly, 
we immediately obtain 

dooi^J-^v) < sup \y - T{y)\ . 
j/esupp(^) 

This comes from ([d]), by using the transference plan t:j- = {\^n x 

Lemma 1. Assume that fi, fl ^ 'P(R^) are two convex combinations: = mono + mi//i and 
fi = moflQ + mifii, where fiQ and Jiq are supported in B{xq, e) for some xq G R^ and e > 0. Then 
doo{fJ',P-) < 2e. 

Proof. Let vri G n(^i,/ii) be the transport plan induced by the identity map, that is 

(/)(x,y)d7ri(x,y) = / (/)(x, x)d/ii(x) 

and let ttq G n(/uo,/io) be any transport plan between /xq and fiQ. Note that vr = movro + mim C 
n(/Li, /i) and supp(7r) = supp(7ro) U supp(7ri). Since vri is supported on the diagonal we have 
sup(3;^y)gsupp(7ri) 1^; — y| = 0. On the other hand, supp(7ro) C supp(^o) x supp(/io) C B{xo,e) x 
B{xQ,e) and therefore sup(^ ,^)gsupp(^|^-) |x — y| < 2e. We conclude that sup(^ j^^gg^pp^^) |x — y| < 2e 
which implies inf^gn(i.,p) sup(^_j^)gsupp(^) |x - y| < 2e. □ 
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3. Lower bound on the Hausdorff dimension of the support 



In this section we consider potentials which are strongly repulsive at the origin and we prove 
that if AW ~ — l/lxl'' as x — >• 0, < /3 < A'^, then the Hausdorff dimension of the support of local 
minimizers of the interaction energy is greater or equal to /?. Actually our result is slightly stronger: 
we prove that if is a local minimizer then the support of any part of fi has Hausdorff dimension 
greater or equal to (3. Let us illustrate the importance of controlling not only the dimension of 
/i, but also the dimension of the parts of /i. Suppose for example that AW ~ — as x — )■ 0, 
then our result implies that any part of fj, has Hausdorff dimension greater or equal to 1. As a 
consequence /i can not have an atomic part. If AW ~ — as x — t- then fj, can not have a 
part concentrated on a curve and so on. 

3.1. Hypotheses and statement of the main result. In this section, we will assume that the 
potential W : — )• (— oo, +oo] satisfies the following hypotheses: 

(HI) W is bounded from below. 

(H2) W is lower semicontinuous (l.s.c). 

(H3) W is uniformly locally integrable: there exists M > such that J^^^ ^-^ W{y)dy < M for all 

X e M^. 

In order to state the main results of this section we will also need the following two definitions: 

Definition 1 (Generalized Laplacian). Suppose W : M.^ — )■ (— oo,+oo] is locally integrable. The 
approximate Laplacian of W is defined by 



-A'W(x) := '^^^^ ( W(x) - / W(x + y)dy 

e V JBiO,e) 



where j-^f^^^ f{x)dx stands for the average of f over the ball of radius r centered at xq, and the 
generalized Laplacian of W is defined by 



-A^W{x) :=liminf|-A(i/")VF(x)|. 



Definition 2 (/3-repulsive potential). Suppose W : — ?■ (— cxo,+oo] is locally integrable. W is 
said to be f3-repulsive at the origin if there exists e > and C > such that 

- A'^Wix) > ^ for all < |x| < e (10) 

- A°VF(0) = +00. (11) 

By doing a Taylor expansion one can easily check that A'^VF(x) = AW{x) wherever W is twice 
differentiable. In particular if W is twice differentiable away from the origin as it is often the case 
for potentials of interest, then (10) simply means that —AW{x) > C/\x\^ for all < |x| < e. 



The terminology "/3-repulsive" is justified by the fact that the rate at which A^W{x) goes to 
— oo as X approaches the origin quantifies the repulsive strength of the potential at the origin, 
therefore the greater /3 is the more repulsive the potential is around the origin. This is the rigorous 
mathematical formulation of what we meant in ([3j). Additionally to hypotheses (H1)-(H3), we will 
need the following technical assumption on the potential W: 

(H4) There exists C* > such that 

A'W(x) <C* Vx G and Ve G (0, 1). 
We are now ready to state the main theorems of this section: 
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Theorem 1. Suppose W satisfies (H1)-(H4) and let fi be a compactly supported local minimizer of 
the interaction energy with respect to the topology induced by doo- IfWis P-repulsive at the origin, 
< P < N , then the Hausdorff dimension of the support of any part of fi is greater or equal to (3. 

Remark 1. Observe that (H3) and (H4) are conditions which restrict the growth of the potential 
and its derivatives at oo. For instance, a potential growing algebraically at oo does not satisfy those 
assumptions. However, if we are only interested in the dimensionality of the support for compactly 
supported local minimizers. Theorem [l] holds under weaker assumptions not restricting the growth 
of the potential at oo. Namely, (H3) and (H4) can be substituted by (H3-loc) and (H4-loc); 
(H3-loc) W is locally integrable. 

(H4-loc) For every compact subset K of there exists Cj^ > such that 

^.'W{x) <Ck Vx G K and Ve G (0, 1), 
with obvious changes in the proof. 

Remark 2. In Theorem^ (resp. Remark^ potential W is assumed to be (3-repulsive at the ori- 
gin and to satisfy hypotheses (H1)-(H4) (resp. (Hl)-(H2)-(H3-loc)-(H4-loc)j. Whereas hypotheses 
(H1)-(H3) (resp. (Hl)-(H2)-(H3-loc) j are easily verified for a given potential, hypotheses (H4) or 
(H4-loc) and the (3-repulsivity are not as transparent. To clarify the meaning of these more technical 
assumptions let us consider the case where W is smooth away from the origin and satisfies 

C 

- AW{x) > -— T for all < < e (12) 



for some < /3 < N . Such a potential satisfies (10) as pointed out in the comment after Definition 
[2| Moreover most potentials of interest satisfying (12) will also satisfy (11) and either (H4) or 



(H4-loc), but of course this need to be checked case by case. In subsection 3.3 we consider some 



typical repulsive- attractive potentials satisfying (12) and we show that they satisfy (11) and either 



(H4) or (H4-loc) depending on their behavior at infinity. 

3.2. Proof of Theorem [l| First note that without loss of generahty we can replace hypothesis 
{HI) by 

(HI') W is nonnegative 

since adding a constant to the potential W does not affect the local minimizers of Ew. The following 
lemma is classical: 

Lemma 2. Suppose W satisfies (HI') and (H2) and let fi G 'P(M^). Then the function : — )• 
[0, +oo] defined by 

V^{x) = {W * /i)(x) = / W{x- y)dfi{y) 

is lower semicontinuous. 

Proof. Suppose then by Fatou's lemma we have 

V^{x) = / W{x - y) dfi{y) < / liminf ^^(xn - y) dn{y) 



< liminf / W{xn - y) d^{y) = liminf F^(x„) 

n J^N n 



as desired. □ 

Suppose now that W satisfies (H1')-(H4). Note that hypothesis (H4) implies that -A°W > -C* 
and as a consequence, for any fi G 'P(M^), the function 

i-A'W * fi){x) = [ {-A'W){x - y)df,{y) = [ [(-AV)(x - y) + C*] dfi{y) - C* 

7 



is defined for all x and -A°PF * : +00]. 

Lemma 3. Suppose that W satisfies (H1')-(H4) and let n e P(M^). // xq is a local min of 
= W * fi, in the sense that there exists eo > such that 

y^ii^o) < ^(2;) for almost every x £ B{xQ,eo), (13) 

then {A"W*fi){xo) > 0. 



Proof. Assume that xq satisfies ( |13| ). We first show that V^(a;o) < +00. If it were not the case we 
would have that = +00 a.e. in B{xq, eo)- But hypothesis (H3) and Fubini's Theorem imply that 



V^{x)dx< f [ W{z)dzdfi{y) < M 



B{xo,l) JM^ JB{xo-y,l) 

and therefore is finite almost everywhere in B{xq, 1), contradicting the fact that = +00 a.e. 
in B{xo, eo). Now, for e < eo we have 



(£ 

/ / W{xo + x-y)dxdfiiy)- [ W{xo - y)dfiiy)] . (14) 
Jr^Jb(o,€) Jr^ I 



< ^i^^±^ I J- v^ixo + x)dx - V^{xo) 



e' 

2(A^ + 2) 

'B{0,e) 

Note that hypothesis (H4) implies that 



j~Bi 



W{xo + X — y)dx < W{xo — y) + 



lB(0,e) " - ^ 2(d + 2) 

Since V^{xo) < +00, the functions y 1— )■ W{xo — y) and y t-^ f^^^ l^(xo + x — y)dx are /U-integrable 
and the difference of the integrals in ( 14 ) is equal to the integral of the difference. Therefore we 
have: 

0< / ^^^^(/ W{xo-y + x)dx-W{xo-y)]dfi{y)= [ A'W{xo-y)dfi{y). {15) 

JrJV e^ \JB{0,e) J JrN 

Because of hypothesis (H4), we have that —A^W + C* > for all e G (0, 1). Therefore using Fatou's 
Lemma and ( 15 ): 

liminf |-A(i/")M/(xo - 2/) + C*| d^(y) 

< liminf / [-A(i/")VF(a;o - y) + C*] d^i{y) < C* , 

that is, {A^W * /i)(xo) > 0. □ 

Proposition 1. Suppose that W satisfies (H1')-(H2)-(H3). Let fx be a local minimizer of the 
interaction energy with respect to the doo o-nd assume that E[^] < +00. Then any point xq G 
supp (fi) is a local minimizer ofV^, in the sense that there exists eo > such that 

y^i^xo) < V^{x) for almost every x G i?(3;o,eo). 

Proof. We argue by contradiction. Assume that there exists xo G supp(/i) which is not a local 
minimum of V^. Fix e > 0. Then there exists a set A C -6(2:0, e) of positive Lebesgue measure, 
such that for x £ A, V^(x) < V^(xo). The set A can be written as follows: 

A = U~ i{x G A; V^{x) < ^^(xo) - 1/n}, 



that is A is an increasing union of measurable sets. Thanks to the continuity from below of the 
Lebesgue measure, it implies that 

< |.4| = lim \{x e A; Vf,{x) < V^{xo) - l/n}|, 

n— >oo 

and there exists uq such that A := {x £ A; < V^{xo) — 1/no} is of positive Lebesgue measure. 

Thanks to the lower semicontinuity of V^, there exists rj S (0, e) such that 

inf V, > V,{xo) - ^ > supF, + (16) 

B{xo,ri) ZriQ ^ ZUq 

Notice that xq G supp(/x) implies fi{B{xo, rj)) > 0. We can therefore define the probability measures 
/^o, by 

MB) = —ji{B n Bixo, v)), l^AiB) = i|S n i| 
mo ^ 1^1 

for any Borel set B E B{MJ^)^ where mo := ji{B{xo, rj)). Let us now write as a convex combination 
/i = mo/^o + Tnim, and define the curve of measures 

fJ-t = {mo - t)fio + tfi^ + mi/xi. 

It is clear by construction that G 7^(M^) for t G [0,mo]. Note that is obtained from /i by 
transporting an amount t of mass from the region B{xo,r]) and by distributing it uniformly in 
the region A. Since both B(xo,rj) and A are contained in B{xo,e), the mass is transported by a 
distance smaller than 2e and therefore we have doo{n, fJ^t) < 2e, see Lemma [T] for details. Inequality 
(16) shows that the function is greater on the region B(xQ,r]) than on the region A, therefore 
one would expect that transporting mass from one region to the other will decrease the interaction 
energy. Indeed we will show that E[fit] < E[fi] for t small enough. Since e was arbitrary, this will 
imply that we can always find a probability measure arbitrarily close to fi (in the sense of the doo ) 
with strictly smaller energy. This is a contradiction concluding the proof. 

We are left to show that Elut] < for t small enough. Since < E[ij,] < oo and given by 

E[fi] = mlElfio] + 2momiB[fio, m] + mjElfii] 

then the three terms E[fiQ], B[fj,Q, fii] and are all positive and finite. Note that -E'[/^^] is also 

finite: indeed, since W is locally integrable by (H3), we have 

■^[/^a] = / / - y) df^li^) df^Aiy) ^ A;^ 1 1 - y) dy < +00. 

JjRNy^^N \A\^ J JB{xo,€)xB{xo,e) 

From (16) and the fact that i3[/z,/xo] < < +oo, we also have that 

B[f^, f^A] + W^< ^[/"' Mo] < +00. (17) 

ZUq 

Using all these, we can show that all combinations of the bilinear form -B[-, •] for the measures no, 
/ii, and fi^ are finite: 

E[iJ,o]<+oo , < +00 , E[fi^]<+oo , B[ni,fio] < +00, (18) 

B[fi^,fio] < —B[fi^,fi] < +00 , B[fi^,fii] < —B[fij^,fi] < +00, (19) 
mo mi 



where we have used (17) in order to obtain (19). Note that in (19) we have assumed mi 7^ 0. If 
mi = then fii can be chosen to be zero and therefore B[fj,^, ni] < +00 trivially holds. Using 



(18)-(19), we are allowed to expand E[fit] as: 
E[fit] = E[{mo - t)fio + mi/ii + tfj.^] 

= (mo - tfE[fio] + mlE[^xi] + t'^E{fi^] 

+ 2(mo - t)miB[fio, Hi] + 2(mo - t)tB[fio, fi^] + 2mitB[fii, fi^] 
= mlE[t,o] + 2mom^B[^,o,^^l]+mlE[^,^] 

+ 2t(moB[mo,fi^] + miB[fii,fj,^]^ - 2t(mQB[Ho, fio] + miB[fiQ, fxi]^ 

+ t^E[^io] + t^E[fi^] - 2t^B[fio, Mi] 

= E[|s] + 2t(B[^l^,^J]-B[f,o,^J])+t^T[|^o,^^A]■ (20) 



Note that in the above computations we have only used the bilinearity of B[-, •] on the space of 
positive measures. However, a formal computation using the bilinearity of -B[-,-] on the space of 
signed measures leads to the same result in a much simpler way: 



E[fit] =E[n-tfio + tfi^] = E[fi] + 2t[B[n^, fi]- B[fio,n]] + t^T[fio, /i^]- 



To conclude the proof note that because of (17) the term B[fi^,fi] — B[fiQ,fi\ appearing in (20) is 
strictly negative and since the term T[fj,Q, /i^] = E[fiQ] — 2B[fio, fi^] + E[nj] is finite we can choose 
t small enough so that E[fit] < E[ij]. This concludes the proof. □ 

Under the additional assumption that W is not singular at the origin, we can obtain a slightly 
stronger version of Proposition [T] which will be needed in section 4. 

Proposition 2. Assume that W and fi satisfy the same hypotheses than in Proposition^ Assume 
moreover that W{0) < +oo. Then any point xq G supp(/i) is a local minimizer ofV^ in the classical 
sense and is constant on any connected compact set K C supp(/x). 

Proof. The proof of the first statement is similar to the proof of Proposition [T] We argue by 
contradiction: assume that S T'(M^) is a local minimizer of E[-] and that there exists xq G supp(//) 
which is not a (classical) local minimum of V^. Fix e > 0, then there exists x^ E B{xo, e) such that 
Vij,{xa) < V^{xq). But since is l.s.c. there exists < rj < e such that 

Vf^ixa) <V^ixo) <Vf,ix) forallxGS(xo,7?). (21) 

We then define /Uq and /ii as in the proof of Proposition [TJ The different idea now is to send mass 
from fiQ to a Dirac Delta at the location x^ instead of distributing it evenly over a set A of nonzero 
Lebesgue measure: instead of letting /ij = (mo — t)fj,o + tfj.^ + mi^i as before, we now define 
fj-t = (mo — t)fiQ + t5xa + rnifii. The same expansion leads to 

E[fit] = E[n-tf,o + t6,J =E[f,] + 2t(B[6x^,H] - B[^lo, l^]) + t^T[fio, 4 J • 



From (21) we obtain that the term — B[fiQ,fi] is strictly negative. In order to conclude 

the argument we need the term T[fj.Q,6xJ = E[fiQ] — 2B[fiQ, 5x^] + -Eif^^a] to be finite. Note that 
i?[(5xt,] = W{0)/2 therefore it is necessary for VF(0) to be finite in order to conclude the proof. 

We now prove the second statement. We follow classical arguments from potential theory, see 
\28\ Proposition 0.4] for instance. Let K he a connected compact set contained in supp(/i) and 
consider the sets A = {x £ K : F^(x) > mix V^} and B = {x £ K : V^{x) = iniK V^}. Since 
is l.s.c. the set A is open relative to K. Let us show that B is also open relative to K. We argue 
by contradiction. Suppose there exists Xa £ B such that for every e > there exists xo,e G K with 
\xa — xo^e\ < e and V^(xa) < V^(xo,e). Then following the exact same steps as in the proof of the 
first statement we can construct a probability measure with lower energy than fi and whose doo 
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distance to /u is smaller than e, therefore leading to a contradiction and proving that B is open 
relative to K. Since K is connected then either A or B must be empty. But since is l.s.c it has 
to reach its minimum on compact sets and therefore A = (/> and B = K. □ 

Remark 3. Since supp(^) is closed, the connected component o/supp(/i) are also closed. So the 
second statement of Proposition [2] implies that Vfj, is constant on any bounded connected compo- 
nent o/supp(/u). In particular if n is compactly supported then is constant on any connected 
component of supp{p). 

Combining Lemma [3] and Proposition [T] we obtain: 

Corolary 1. Suppose that W satisfies (H1')-(H4). If fi is local minimizer of the interaction energy 
with respect to d^o and E[fi\ < +oo, then (A*^!^ * ^ for all x £ supp(/i). 

We recall the following result from [16, Theorem 4.13]. 

Proposition 3. Let A be a Borel subset ofMJ^, and s > 0. // there exists a probability measure 
^ G ■p(M^) supported on A such that 

d^{x)d^{y) 

— I 1 < 

then dimnA > s, with dimn being the Hausdorff dimension of A. 

We are now ready to prove the main theorem. 
Proof of Theorem^ Let p be a nonzero part of /x, that \s = p + u for some nonnegative measure 



V. Let A = supp(/9) and let us show that dminA > /3. Choose e small enough so that (10) holds, 
choose xq £ A and define the measure 

fio{B) = p{BnB{xo,e/2)). 

Clearly p can be written p = pQ + pi, where po and pi are two (nonnegative) measures of mass 



mo > and mi > and where po is supported in ^4 n B{xq, e/2). Then from (10) we get: 



C 1 1 ^J^p^ < / / _aOt^(x - y)dpoix)dp,iy) 



\x - y\ 



XJl 



-^'^W{x-y)dp{x)dpQ{y)- // -/\''W{x-y)dpi{x)dpQ{y) 
(A^W^ * p){y)dpa{y) + [[ A^W{x - y)dpi{x)dpo{y) 

< // A°W{x - y)dpi{x)dpo{y) < C*mimo < +00. 

We have used the fact that A^VF * ^ is nonnegative on the support of p from Corollary [l] and that 
A^W{x) < C* by hypothesis (H4). We then apply Proposition [s] to the probability measure po/mo, 
which is supported on A, to obtain dimnA > f3. □ 

3.3. Example of potentials satisfying the hypotheses of Theorem [ij In this subsection we 
consider the class of potentials 

Wa{x) = Cha{x) + Tpix) 

where ip G C'^(M^) bounded from below, c > and ha '■ — )• (— oo, oo] is the power-law function: 

ha{x) = —\x\"/a 

for x ^ and a € M with the convention ho{x) = — log|x|. We define ha{0) = if a > or 
ha{0) = +00 if a < 0. The potentials Wa are typical examples of repulsive-attractive potentials 
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behaving like — |x|"/a around the origin. It is trivial to check that Wq. satisfies (Hl)-(H2)-(H3-loc) 
for any a > — (in the case a > the function tp need to grow fast enough at infinity for hypothesis 
(HI) to hold). Note also that for x ^ we have 

- AW^{x) = ^ ^ - A^{x) (22) 



and therefore ifa + A^ — 2 > then Wa satisfies ( 10 ) from the definition of /3-repulsivity with 
P = 2 — a. The goal of this subsection is to show that Wa also satisfies (11) and (H4-loc). 
We start by checking (|11[). An explicit computation gives 



V Tbim I \+oo if2-Af<a<0 

where we have used the fact that /iq,(0) = for a > and /iq,(0) = +cxd for a < 0. Letting e — )• 
and using the fact that A'0(O) is finite we obtain 

- A°VF(0) = +00 for all a < 2. (23) 



Combining ( |22[ ) and (23) we see that for 2 — < a < 2 the potential Wa is /3-repulsive with 
/3 = 2-aG (0,iV). 

We now show that for a > 2 — N the potentials Wa satisfies hypothesis (H4-loc). The key 
point is that the functions ha are superharmonic for a > 2 — N. Let us recall the definition of 
superharmonicity: 

Definition 3. A lower semicontinuous function h : — t- (—00,00] is said to be superharmonic 
on the connected open set Q if it is not identically equal to +00 on $7 and if 

h{x) > I h{y)dy 

JB{x,r) 

for all X and r > such that B{x, r) C il. 

We also recall that if h £ C^(r2), then h is superharmonic on O if and only if Ah{x) < for all 
X £ Q. To see that the functions ha are superharmonic for a > 2 — N , first note that for x 7^ 

. , .X (a + A^-2) 
A/i„(x) = -^^^r^<0- 

Therefore ha is superharmonic on M^\{0} and it can be easily checked that it satisfies the super- 
mean value property at the origin Definition 2.1]. Both together imply that it is actually 
superharmonic on the full space M.^ , [28', Corollary 2.1]. As a consequence we directly obtain from 
the definition of the approximate Laplacian that A'^haix) < and therefore 

A'Wa = cA'ha + A'V < A'V'- 

To conclude we note that since ip G C^(M^) a simple Taylor expansion shows that A'^ip converges 
uniformly to Aip on compact sets. Indeed, the expansion gives 

2(N + 2) f 

A'tl;{x) = ^^ — -f Tpix + y) - ^{x) dy 

e JB{0,e) 

2(N + 2) 

^ ^'1 y'V^{x)+y''H^{x)y + 0{e'')dy (24) 




yfdy + 0{e^)] (25) 

B(0,e) J 

Ai;{x) + 0(e) 
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To go from (24) to (25) we have used the fact that most of the terms in the Taylor expansion are 
equal to zero after integrating on spheres due to symmetry. The only remaining terms are the ones 
involved in the Laplacian. Note that since the partial derivative of order 3 of V' St.re bounded on 
compact subsets of M'^, then the error term is uniform for x in compact sets. Since A^/^ is bounded 
on compact sets, we conclude that for a > 2 — the potential Wa satisfies (H4-loc). 

Moreover, let us point out that if ip is well behaved at infinity in terms of regularity and growth, 
then Wa satisfies either (H1)-(H4) or (Hl)-(H2)-(H3-loc)-(H4-loc). We summarize this discussion 
in the following proposition: 

Proposition 4. If2 — N<a<2 and if ip £ C^(]R^) then Wa is (2 — a)-repulsive around the 
origin and satisfies (H4-loc). 

4. Mild Repulsion implies 0-dimensionality 

In this section, we will show that if the potential is mildly repulsive, meaning that it behaves 
locally near zero like — with a > 2, then the support of the measure cannot contain measures 
concentrated on smooth manifolds of any dimension except 0-dimensional sets. 

Definition 4. Let 1 < k < N . A probability measure n G T'(M^) is said to have a regular 
k- dimensional part if it can be written 

/i = /il + //2 

where /ii is a nonnegative measure on and fi2 is defined by 

i;{x)dfi2{x) = [ ^{x)f{x)da{x) V?/; G Co(M^) 
JM 

for some manifold Ai of dimension k and a non identically zero nonnegative function / : — )• 
(0,+oo] integrable with respect to the surface measure da{x) on Ai. Moreover to avoid pathological 
situations, we assume that there exists xq €z A4, c, k > such that 

f{x)>c yx e MnB{xo,K). (26) 

We now state the main result of this section: 

Theorem 2. Let W G C7^(M^) be a radially symmetric potential which is equal to — jxl'^/a in a 
neighborhood of the origin. If a > 2 then a local minimizer of the interaction energy with respect 
to doo cannot have a k-dimensional component for any 1 < k < N . 

For the above theorem to be true it is not necessary for the potential to be exactly equal to a 
power law — Ixl*^, a > 2, around the origin. It is enough for the potential to behaves like — 
a > 2, at the origin in a precise convexity sense - see Theorem [3| 

4.1. Preliminaries on convexity. To prove Theorem [2| we need some convex analysis concepts, 
see \13\ [1] and the references therein. The term modulus of convexity refers to any convex function 
on the positive reals satisfying 

{(po) (j) : [0, oo) — > M is continuous and vanishes only at 0(0) = 0. 

{(f)i) (j){x) > —kx for some k < co. 

Now, we can quantify the convexity of certain functions in terms of a modulus of convexity. 

Definition 5. A function h : [0, +oo) — )■ M is cp— uniformly convex on (a, b) if there exists a modulus 
of convexity (j) such that 

h ("^^V^) < I {h{n) + Hr2)) - ] <A(t) dt, (27) 



2 J - 2' ' ' ' " 4 

for all ri,r2 G (a, b). 
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A function h : [0, +cxd) — )• M is \— convex on (a, b) if it is (j)-uniformly convex with (j){s) = \s and 

A e M. 



Note that if h is A— convex, then (27) reads 



h I^^Y^J ^ 2^(^i) + 2^(^2) - 8 - r■2)^ (28) 

for all ri,r2 G {(i,b). It is equivalent to assume that the function h(r) — ^r'^ is convex on {a,b). 
The following proposition can be easily proven: 

Proposition 5 (Convexity properties of power laws). 

(i) If q e (1,2], then h{r) = is X- convex on [0,R] for A = inf(o,_R) h" = q{q - l)i?'?"2 > 0, 

and thus, uniformly convex on [0,R]. 
(a) If q > 2, then h{r) = r'^ is (p— uniformly convex on M4., with (j)(t) = 'i^^'^t'^^^ j q. That is 

h (^^4^) - I (^('■1) + ^('■2)) - 2-''\n - rsr, (29) 

for ri, r2 > 0. 

4.2. Proof of Theorem [2| In this subsection we prove the following generalization of Theorem[2j 

Theorem 3. LetW{x) = w{\x\) be continuously difjerentiable, bounded from below, and decreasing 
as a function of\x\ in a neighborhood of the origin. Assume moreover that W behaves like the power 
law — a > 2, near the origin, in the sense that for some C* > and R > small enough, 
W{r) = —h{r'^) satisfies: 

• ifae (2,4], h is X-convex on [0,R] with X = C*i?°/2-2_ 

• if a £ (4,00), h is (p— uniformly convex on [0,R], with (f){t) = C*t"/^^^, 

and C*\w'{r)\ < r°^^^ on [0,R]. Then a local minimizer of the interaction energy with respect to 
doo cannot have a k-dimensional part for any 1 < k < N . 

Theorem [2] is a direct consequence of Theorem |3j thanks to Proposition [5| 

We first provide an explicit formula for how the energy changes when perturbing a local mini- 
mizer: 

Lemma 4. Suppose that W : — )■ (— cxd,+oo] is symmetric, l.s.c. and bounded from below with 
W{0) < +00. Let jj G 'P(M^) be a local minimizer of the interaction energy with respect to doo 
and E[i^i] < +00. Given a connected domain Q C supp(/x), a Borel map tt : — )■ and a convex 
decomposition fj, = mi/xi + m2jU2 with supp(^i) C il, we deduce that 

E[mi{TT#fii) + m.2H2] - E[miiJ.i + 711,2/^2] = T[it#hi,ij.i] . (30) 

Proof. Since in ([T]) is a bilinear form, 

E[mi{TT#^ii) + m2H2] - E[mifii + 7712^2] 

= mlB[TT#fii,TT#ni] + 2mim,2B[TT#ni,fj,2] 

-mlB[fj.i,fj,i]-2mim2B[fj.i,fi2]. (31) 

We now use the fact that is local minimizer to express the terms involving /i2 as terms involving 
only fii. Proposition [2] implies that the function V^{x) is constant on the connected domain and 
since 7r(r2) C 1^ we have: 



WiTr{x)-y)dfi{y)= [ W{x-y)dfi{y) Vx G Jl. 

Jrn 
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and therefore, since /u = mifii + m2fJ.2, 



mi / W{Tr{x) - y)dfii{y) +7712 / W{Tr{x) - y) dfi2iy) 



= mi W{x - y) dfii{y) + m2 / W {x - y) dfj,2iy) 

for all X G fi. Since supp(//i) C 0, we can integrate both sides with respect to d^i{x) and obtain, 
after multiplication by mi: 

m\ I W{tt{x) -y)dp,i{y)dfii{x) + 17111712 W {'i:{x) - y) diX2{y) d^i{x) 

= 7n\ I W{x-y)dfj,i{y)dfj,i{x)+mim2 W{x - y) dfj.2{y) dfj.i{x) 
JM^xR^ Jr^xR'V 

or equivalently, using the i?-notation, 

2m^i?[^i, 7r#/xi] + 2mim2i?[;U2, vr#^i] = 2mf i3[^i, ^i] + 2mim2-B[^2, ^i] 

and therefore rearranging the terms, we can express the terms involving fi2 in terms of the ones 
involving only /ii: 



2mim2 



B[n#ni,fi2] - B[ni,fi2] 



2mt 



The desired identity (30) is readily obtained by plugging the last equality into (31) reminding the 
definition of T[/u, i'] in Section 2. □ 

Definition 6. Let 1 < k < N . We denote by the the k-dimensional disk of radius e: 



Dt 



'e — {(^1; ■ ■ • > ^A^) • xl + . . . + xl< and Xk+i 
and by Ve,k £ 'P(M^) the uniform probability distribution on : 

1 ^ 



XN = 0}, 



where \D^\ is the Lebesgue measure of dimension k of D^, that is \D^\ = a^e^ where the area 
of the unit k-dimensional ball, v^^^ then satisfies 



il){x)dv^^k{x) 



1 



il)(xi, . . . ,Xk,Q, . . . ,Q) dxi . . . dxk 



for alliPe C°(M^). 



The following Lemma combined with Lemma|4]shows that if a flat A:-dimensional disk is contained 
in the support of a local minimizer, then the energy can be reduced by concentrating all the mass 
contained in the disk into a single point. As a consequence the support of a local minimizer cannot 
contain a flat fc-dimensional disk. 

Lemma 5. Suppose that W{x) = — satisfies the assumptions of Theorem^ Then, there 
exists Ck^a > such that for s small enough, 

T[6o, Ve^k] = B[6o, 5o] - 2B[5q, Ve,k] + B[Ue,k, l^e,k] < -Ck,ae". 
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Proof. Since W is bounded from below and W{0) < +00, we can assume without loss of generality 
that W{0) = by adding to a suitable constant. Then, the first term B[6o, do] is equal to zero. 
Symmetrizing the integral involved in second term we obtain: 



K\y?)dve,k{y) 



h{\y\ ) dve,k{x)due,k{y) 



Since Vg^ki—y) = '^e,kiy), we can also symmetrize the third term and obtain: 

1 
2 
1 

/R^xRJV 



B[l^e,k,l^e,k] 



h{\x - yp) dve,k{x)dve,k{y) 
[h{\x - + h{\x + dvg^k{x)dve,k{y)- 



Combining the three terms we find 



T[5o,Ve,k\ 



1 



R^xRJV 



A{x,y) dve,k{x)dve,k{y) 



(32) 



with 



A{x,y) ■.= h{\x\') + h{\y\') 



2. h{\x - y\'^) + h{\x + y\ 



Under the assumptions of Theorem [2| /i is convex on (0,2e2) and since h{0) = 0, we deduce 

^2 



Hrf) < 



for > 0, i = 1, 2. Using the above inequalities for i = 1, 2 we get 

/I 1 

h{\x\^) + h{\y\^) < h{\x\^ + = h{-\x + y\^ + -\x-y 

In the rest of this Lemma, C will denote some generic constant that will change from step to 
step. For a G (2,4], h is A— convex on (0, 2e^) with A = Ce"~^, so that plugging into (28), we 
obtain 



/i(|x|2) + /i(|y|2) < \ (|x + + \ {\x - - Ce'^-\x ■ yf. 



(33) 



Combining (32) and (|33j) we get: 

T[5o,v^^k\<-Ce°'~'^ I {x-yf dve,k{^)dve,k{y) 

JR^xR^ 

= -Ce" / (x • yf dui^k{x)dui^k{y) ■ 

JR^xR^ 

For a > 4, /i is convex with (j){t) = Ct"/^~^, so that plugging into (|29[), we obtain 



/i(|x|2) + /i(|y|2) < ^/i {\x + y|2) + ^/i (|x - - C|x • 2/1-/2. 



(34) 



Combining (32) and ( |34[ ) we get: 

T[6o,iy£,k] < -C \x-y\°-/'^ dve^k{x)due,k{y) 

JR^xR^ 



RJ^xR^ 



• yl"''^ dvi^k{x)dui^k{y) ■ 



□ 
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The last Lemma combined to Lemma |4] shows that the support of a local minimizer cannot 
contain a flat A;-dimensional disk of radius e. To conclude the proof of Theorem [3j we need to 
introduce some differential geometry tools. Let i? > 0, and g : — )• a C^-function such 

that g{0) = 0, Vg{0) = 0. We define the parameterisation Pg of the graph of g as follows: 

Pg-.Dii M^, (35) 

(x',0) H> {x',g{x')) 

where x' = (xi, . . . ,Xk) stands for the k first coordinates. Let us remark that classical differential 
geometry implies that any C^-manifold can be locally parameterized by such graphs by choosing 
conveniently the axis and reordering of variables. Moreover, this can be done in such a way that 
the volume element of the graph Jg is as close to the unit volume element of the fiat tangent space 
by taking R small enough, see [22] • More precisely, there exists a constant Cg depending only on 
the second derivatives of g on such that 

\\Jg - M\l°°{DI) ^ CgE, (36) 

for < e < R small enough. 

Lemma 6. IfW satisfies the assumptions of Theorem^ and g G (7^(1?^, M^"*') satisfies g{0) = 0, 
V5(0) = 0, then for e > small enough, 

Proof. Note that by continuity for e > small enough, we have HV^H^^oo < 1. We first point 
out that 



with 



T[6o,Pg#l^e,k]-'^[^0^'^e,k]= / Mx,y)dl^e,ki^)dl'e,kiy) 

, w(\P„(x) — P„(y)\) — w(\x — y\) r ,, „ , 
A{x,y) = ^' ^ ^ ^- w{\Pg{x)\) - w{\x\) 



2 

Thanks to the definition of the parameterisation Pg, \Pg{x) — Pg{y)\ > |x — y|. Moreover since w 
is decreasing in a neighborhood of the 0, the first term in A{x,y) is negative for max(|x|, \y\) < e 
small enough. To estimate the second term, we use the mean value theorem for g around x' = 0, 
remembering that g{0) = 0: 

\Pg{x) - {x',0)\ = \g{x') - g(0)| < e\\Vgh^^o,), 
since C*|w'(r)| < r^~^^ we conclude 

w{\Pgix)\) - W{\x\) < |h'||Loo([o,2£])e||V5r||ioo(£,fe) < II Vff ||z.cx.(^fc) . 

□ 

Proof of the Theorem^ Assume that ;U is a local minimizer of E in doo and that it has a regular 
/c-dimensional part in the sense of Definition |4] Let M. be the C^-submanifold on which this 



component is supported, and / be the density on M of this component. Due to assumption (26) 
there exists xq G , c, k > such that 

f{x)>c, yx £ MnB{xo,K). 

As discussed above and without loss of generality, we can assume that xq=0 and that M is locally 
the graph of a C^-function g : D'^ — > M^~'^, for some k> R> 0, such that g{0) = 0, V5'(0) = 0. 
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Let Pg be the parameterisation defined in (35). Note that for e < R, fif := Pg^u^^k £ P{ 



IS 

absolutely continuous with respect to the volume element on Ai with a density denoted still by /if 
satisfying 

M\\L^iM4^) < ^ < \Dl\^i_Cge) ' "^'^ = Ml '^^^^^^ 

where we used (36). Therefore, choosing mi = ^\D^\ (1 — Cgs), then f{x) > mi/if on x G D^, and 
we can decompose /x as a convex combination 



mini + "^-2/^2) 



where |U| G ■p(M^). 

We are going to send now all mass from fi^ to a Dirac Delta at xq = 0. Let us define tt : — > 
by vr = and fi^ := mi7r#/if + m2H2, then vr#/i| = Jq- Moreover, /i^ is a small perturbation 
of H in (ioo- 

dooif^,l^')<e{l + \\S/g\\Loo^ni))- (37) 

To check this just take a map 7" in Definition [o] such that T{x) = x for all x G M./Pg{D^) and 
such that r(x) = for x G Pg{D^). Thus, the maximum displacement produced by the transport 
map T is bounded by the maximum of \Pg{x)\ for x G P3(L>^) leading to ([37]) using that g(0) = 
and the mean value theorem. 

Since ii\ has a connected support that contains 7r(supp(/Ltf )) = {0}, we can apply Lemma |4] to 
get: 

E[^^']-E[^l] = mlT['IT#^,\,^l\] 

= ml T[7T#fil, + ml (T[7r#/if , /.f] - r[^#/.f , v,^k]) 

Since vr^^/xf = 5q, we can use Lemma[5]to estimate the first term, and since moreover /if = Pgij^i'e,ki 
we can use Lemma [6] to estimate the last two terms, so that we finally conclude 



E[ix'] - E[,A < ml -Ck,^e^ + Ce"||V<7k. 



'(Dl) 



Since g G C-'^(-D^) and Vg{Q) = imply that || V(7||j;^cx)(£)ft) — > as e — ^ 0, thus if e > is small 
enough, E[pf \ - S[/i] < 0. 

Thus, /i^ is a better competitor in the minimization of E for e arbitrary small. This leads to a 
contradiction with the fact that /i is a local minimizer of E showing Theorem [2j □ 



5. EuLER-LAGRANGE approach to study local MINIMIZERS in the (i2-TOPOLOGY 

So far we have used transport plans to build perturbed measures. This enabled us to study local 
minimizers of the interaction energy with respect to the doo-topology. To study local minimizers 
with respect to the (i2-topology it is actually possible to use a more classical Euler-Lagrange ap- 
proach as we will present in this section. The Euler-Lagrange conditions that we will derive were 
formally obtained in [3] by perturbing densities inside and outside their support. Here, we provide 
a fully rigorous proof in the case of probability measures endowed with the distance d2- 

Theorem 4. Given an interaction potential W satisfying (H1)-(H2). Let us consider // G V2{^^) 
a local minimizer of E with respect to d2 such that E[ij] < oo. Then, 

(i) {W * ^x){x) = 2E[^] fi-a.e. 
(a) {W * /i)(x) < 2E[fi] for all x G supp(/i). 
(Hi) {W * /i)(x) > 2E[fj] for a.e. x G M^. 
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Proof. As usual, we assume that > without loss of generality. Lemma [2] implies that W * is 
well defined, lower semicontinuous, and non-negative. 

In order to prove the first two items, let us choose (p G C^(R^) to define 

99 — / ipd^ I /u := a(x)/i 
Jrn ) 

and [i^ = [i^tv = (1 + ea(x))/x with e > to be specified. It is clear that ;Ue(M^) = 1 since a{x) has 
zero integral with respect to [i. Moreover, since o(x) > — 2||93||icx) then /x^ > for e < ^y^y^^^ = e^. 

Therefore, /Xg e P(M^) for all e < e,^. It is easy to check that ii^ G ^2(1^^)1 that fi^ ^ fi weakly-* 
as measures, and 



In fact, since ea(x) converges uniformly to 0, these claims follow by dominated convergence theorem. 
Therefore, we conclude that 

(i2(/ie, /i) — 5- as e — 0. 

Note that it is not true that doo{f^e,f^) — as e — t- since for localized test functions (p in subsets 
of the supp(//), we are always forced to move mass in supp(/x) for a fixed distance not depending 
on e. 

Now, since /i is a local minimizer in d2 then E[fj,f] > E[ij] for e small enough. Moreover, since fi 
has finite energy, then E[fj,f] < 00 and we can expand it as 

^^^^^"^^^^ = // W{x- y)du{x)dfi{y) + [ [ Wix - y)dy{x)dv{y) > , 

with both integral terms well-defined. As e — )• 0, we easily get 

W{x - y)dv{x)d^x{y) > 



or equivalently, 

V9 [{W * n){x) - 2E[fi]] dfi{x) > 



for all if S Cq°(M^). Since one can take either cp or —cp as test functions, we deduce 

j ip[{W* fi){x) - 2E[fi]] dfiix) = 

for ah if e C^{R^), and thus (i) is satisfied a.e. /x. 

Let us now prove (ii). Take x £ supp(^) then there exists {xn}neN — ^ x with x„ G supp(/i), such 
that {W * ^j){xn) = 2E[^]. The existence of such a sequence is ensured since fi{B{x,e)) > for all 
e > by definition of the support of /i. Then, by lower semicontinuity of W * /j, we get 

{W * ii){x) < liminf(M^ * n){xn) = 2E[fi] . 

and then (ii) is satisfied. 

In order to show (iii), we consider different variations to the ones constructed above. Take 
tp G C^{R^), > and then take 

v = Tp — I / tp dx] fl. 

Again, defining fif, = fi + ev, then it verifies fi^i^^) = 1 and if e < 1/ J tp dx then /x^ > 0. Let us 
remark that this cannot be done for a changing sign test function ijj. As previously, it is easy to 
check that 

(i2(/ie. A*) — >■ as e — ^ 0, 
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note that it is not true that dooipe, fJ-) — )• as e — )• since we always need to transport some mass 
from outside the support of fi to supp(/i). 
Proceeding similarly as in point (i), we get 

W{x - y)du{y)dfi{x) > 

taking e — )■ in E[ij,^] > E[fj]. Therefore, we conclude that 

{{W * fi){x) -2E[fi])tl>dx > 0, 



for ah V G C^(M^), i>>0. This readily impHes (iii). □ 
Remark 4. Note that putting together (i), (ii), and (iii) in previous theorem, we conclude that 

{{W * = 2E[^] for a.e. x G supp(//) 

{W*fi){x) >2E[fj] for a.e. x G \ supp(//). 

if is absolutely continuous with respect to the Lebesgue measure. These two properties are the 
Euler- Lagrange conditions that were found for densities in [4j . 

Remark 5. Let us now clarify the differences between local minimizers in the d2 — topology and 
local minimizers in the doo — topology. Following |17|, let consider as an example the interaction 
potential W{x) := —x + ^ in one dimension. Then, 

Pm = m6o + (1 - m)6i 

is a critical point of the interaction energy for any m G [0, 1]. Theorem 3.1 in [T7j shows that the 
measure pm is a local minimizer in the doo — topology o,s soon as m £ (1/3,2/3). Lndeed, what is 
proven is the stronger statement that pm is locally asymptotically stable for the aggregation equation 
Q with respect to any perturbation in the doo — topology. However, E{pm) = 5 ("^ ~ 5) ~ 
that px/2 only can be a local minimizer of the energy in the d2 — topology (and one can prove that it 
actually is). This shows that the set of local minimizers with respect to the d2— topology is strictly 
contained in the set of local minimizers with respect to the doo — topology. Moreover, numerical 
simulations suggest that, form G (1/3,2/3), pm is actually stable (although not asymptotically 
stable) with respect to small d2— perturbations. As a consequence, when using a gradient flow 
approach to compute numerically minimizers of the interaction energy via particles, one obtains 
doo — local minimizers which typically are not d2 — local minimizers (see e.g. Fig 2 of |17j ). 

6. Numerical Experiments 

In this section we conduct a numerical investigation of the local minimizers of the discrete 
interaction energy ([2| with high number of particles. The gradient flow of ^ is given by the 
system of ODEs: 

n 

Xi = -Y^ mjVW{Xi - Xj). (38) 
j=i 

In order to efficiently find local minimizers of ([2]), we solve (38) by an explicit Euler scheme with 
an adaptive time step chosen as the largest possible such that the discrete energy Q decreases. 
This scheme is nothing else than a gradient descent method for the discrete energy Q. Although 
this method might not be accurate enough for the dynamics, it is efficient to find local minimizers 
of the discrete energy. In stiffer situations an explicit Runge-Kutta method is used instead. These 
methods are essentially the ones proposed in [361 ES] . 
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Table 2. Minimizers of Ey^ in 



for various power-law potentials with n = 2, 500. 



The results of these simulations in two dimensions with power-law potentials were presented in 
the introduction, see Table [T} In Subsection we discuss similar computations in three dimen- 
sions. We also provide numerical experiments suggesting that for some potentials, there are local 
minimizers of the interaction energy with mixed dimensionality, that is, local minimizers that are 
the sum of measures whose support have different Hausdorff dimension. 



In Subsection 6.2 we show how our numerical results can be further understood by using the 
results from ^23 ^ 136 ^ 13]. where a careful stability analysis of ring solution (in 2D) and spherical shell 
solution (in 3D) was conducted. We also show how this stability analysis connects to the analytical 
results presented in this paper. 



6.1. Numerical experiments in 3D. We first compute numerically local minimizers of Ey^ where 
W is the power-law potential defined by ([s]). Recall that lS.W{x) ~ — l/|x|^ with /3 = 2 — a as 
X — 7- 0. The computations are performed with n = 2, 500 particles. The results are shown in Table [2] 
and are discussed below: 
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• Subfigure (a): a = 2.5 and 7 = 5. The support of the minimizer has zero Hausdorff 
dimension in agreement with Theorem [2} Actually, in this particular case it is supported 
on 4 points forming a tetrahedron. 

• Subfigure (b) and (c): the two potentials have the same behavior at the origin, a = 1.25, but 
different attractive long range behavior (7 = 15 and 7 = 1.4 respectively). Theorem [T] shows 
that the Hausdorff dimension of the support must be greater or equal to /3 = 2 — a = 0.75. 
Numerically, we observe that the local minimizer for the first example has a two-dimensional 
support and the minimizer for the second example has a three-dimensional support. We did 
not choose the value a = 1.5 because we were not able to obtain a change of dimensionality 
of the stable steady states varying 7 > a. Note that a = 1.5 is always above the instability 
curve for radial perturbations which meets line a = 7 at the point (\/2, \/2)- See Figure 
and Subsection 6.2 for more details. 

• Subfigure (d) and (e): the two potentials have the same behavior at the origin, a = 0.5, but 
different attractive long range behavior (7 = 23 and 7 = 1.4 respectively). Theorem [T] shows 
that the Hausdorff dimension of the support must be greater or equal to/3 = 2 — a = 1.5. 
Numerically, we observe that the local minimizer for the first example has a two-dimensional 
support and the local minimizer for the second example has a three-dimensional support. 



2000 4000 6000 8000 10000 

Particle Numbor 



Figure 1. Distances of the particles from the center of mass for the power law 
potential with a = —0.5, 7 = 5 in 3D. Case (f) in Subsection 6.1 in Table [2] with 
n = 10,000. 

• Subfigure (f): a = —0.5 and 7 = 5. Theorem [T] proves that the Hausdorff dimension of the 
support must be greater than P = 2 — a = 2.5, which can also be observed numerically. In 
Figure [T| we have represented the radius of particles to the center of mass. The particles 
seem to organize into successive two dimensional layers. Such lattices were also observed in 
[22], and it is related to the finite number of particles used in the simulations. 

Notice that we were not able to find examples of interaction potentials leading numerically to a 
local minimizer with one dimensional support. We could however observe such situations with an 
additional asymmetric confining potentials, we thus believe it should be possible to produce such 
cases. 

A natural question following Tables [l] and [2] is whether it is possible to produce local minimizers 
that are a sum of two measures whose support have different Hausdorff dimensions. A possible 
candidate was already observed in [36]. Here, we analyze it more carefully with much larger number 
of particles. From our simulations, it seems that the interaction potential W{x) = wdxl) with w 
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defined by 

—w'{r)=tanh.{{l — r)a) + b, a = 5, 6 = 0.5, 

leads numerically to a local minimizer consisting in a ball (Hausdorff dimension three) inside a 
spherical shell (Hausdorff dimension two), see Figure [2j 





2000 4000 6000 

Particle number 



Figure 2. Left: Local minimizer in 3D with n 
particles from the center of mass. 



10,000. Right: Distance of the 



The distance of each particle to the center of mass is displayed on the right part of Figure [2j 
The inner ball appears to be composed of five equally spaced layers of particles. This is most likely 
due to the fact that particles are organized into a lattice configuration, and therefore the distances 
between the particles and the origin do not form a continuum. It is instructive to compare the 
distribution of the radius of the particles in the right subplot of Figure [2] with the one in Figure [T] 
for the case of an approximated local minimizer with three dimensional support, i.e.. Case (f) of 
Table [2] Although Theorem [T] guarantees that the support of the local minimizer corresponding to 
Figure [T] has Hausdorff dimension greater or equal to 2.5, we can also observe that particles arrange 
themselves in layers. Notice that in dimension N=2, such artifacts also appear in simulations using 
a finite number of particles, see Figure 4 in [22j . 

6.2. Relationship with previous works on ring and shell solutions. An important charac- 
teristic of the analysis performed in the main theorems of this work is that we do not assume a 
specific shape on the local minimizers. If on the contrary, one is interested by the special case of 
delta ring minimizers (in 2D), or spherical shell minimizers (in 3D), perturbative methods provide 
more detailed results. 

In [23] the local stability of discrete ring solutions, made of A^-particle equally distributed in a 



circle, was studied for the A^-particle system (38). The authors considered the power law interaction 
potentials ([s]), and led a formal linear stability analysis for the continuum ring solution as steady 
state of the aggregation equation ([g]) by taking N ^ oo. Those predictions were then confirmed 
numerically. They could not obtain nonlinear stability of the ring solution particularly because there 
is no spectral gap as — t- oo, i.e, the largest negative eigenvalue tends to when N — t- c«. In 0, 
the nonlinear stability of the ring solutions was proved for radial perturbations, corroborating some 
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of the formal results of |23] . together with the instability due to fattening in the complementary 
set of parameters. 
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Figure 3. Sketch of all the computed cases in dimension N = 2. The parameters 
(7, a) used in Table [l] are marked with *, while those used in Table [s] are marked 
with o. Notice that q < 7 is necessary for the interaction potential to be confining. 
In dark gray is represented the set of parameters such that a delta ring could be a 
local minimizer. 

We have represented this set of parameters in Figure [3j as well as all the parameters used in the 
two dimensional numerical simulations of this article (Tables [T] and [s]) . As the caricature presented 
in Table [3] shows, crossing the lower border of this set, curve a = 7/(7 — 1), leads to a "fattening" 
of the delta ring, that is to minimizers with dimensionality 2, see |23l |3]. On the other hand, 
crossing its upper border, given by the curve marked with 3, does not modify the dimensionality 
of the stable steady states as long as a < 2 (they remain one dimensional), but leads to a "shape" 
instability towards a triangular configuration that breaks the ring into 3 connected one dimensional 
components as in case (b) of Table [ij 

Finally, if a > 2, local minimizers become of dimensionality 0, as predicted by Theorem [2| 
whereas if a < 1, all the minimizers are of dimensionality 2, as shown by Theorem [T] 













a = 0.01 


a = 1.1 


a = 1.5 


a = 2.2 



Table 3. Evolution of local minimizers when a > increases, while 7 = 5 remains 
constant. The computations were done with n = 10, 000 particles. 

In three dimensions, a linear stability analysis of discrete spherical shell solutions is also possible 
but it leads to more cumbersome instability curves, see [361 ES]- Again, the results in |3] give 
the "fattening" instability curve dividing instability from stability under radial perturbations. In 
Figure [4j we have only represented the set of parameters such that the spherical shells are not local 
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minimizers for spherically symmetric perturbations, as well as all the parameters (7, a) used for 3D 
numerical simulations in this article in Table [2j Just as we have observed in the 2D case, crossing 
the lower border of this set leads to a "fattening" instability of the spherical shell. 




Figure 4. Sketch of all the computed cases in dimension = 3. The parameters 
(7, a) used in Table [2] are marked with *. Notice that q < 7 is necessary for the 
interaction potential to be confining. The curve is the limit between parameters 
leading to spherical shell solutions for radial perturbations (above the curve) and to 
minimizers of dimensionality 3 (below the curve) . 



Notice finally that it is also possible to modify the dimensionality of the local minimizers with 
other perturbations of power law potentials. As an example, in Table [4j we consider the following 
perturbations of the power law potential ^ : 

U|° Ixl''' 3 

W(x) = 1 1 cos(px), a < 7, p = 3,5. (39) 

a 7 2p 

In Table [4] we have represented the power-law case in the first column, and the perturbations 
in the next two. For (7,0;) = (2, 1.5), the unperturbed power-law potential leads to a local mini- 
mizer with Hausdorff dimension two. When we add the perturbation p = 3, the dimension of the 
minimizer changes to one. Notice that the perturbation does not alter the local behavior of the 
potential at the origin or at infinity, suggesting that Theorem [T] is probably sharp at least in terms 
of natural dimensions. 
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Table 4. Local minimizers with the power-law potential ([s]) and the perturbed 
potential (pi), n = 10,000. 
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